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ABSTRACT 


The  moaning  of  and  need  for  a  "unipotential"  for  the  valence  elec¬ 
trons  of  an  atom  are  presented,  with  a  treatment  of  an  ad  hoc  proposition 
due  to  J.  C*  Slater  for  such  a  potential* 

2 

The  Hartree-Foclc  field  for  the  configuration  Is  2s2p:?S  of  carbon 
is  evaluated  from  radial  distribution  amplitudes  given  by  A*  Jucys,  and 
tabulated  in  such  a  manner  as  to  facilitate  unipotential  calculations* 

The  first  phase  of  a  self-consistent  calculation  on  Slater's  unipotential 
is  then  effected* 

The  accuracy  and  immediate  applicability  of  the  resulting  radial 
distribution  amplitudes  are  then  discussed* 
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Glossary  of  Symbols  not  Defined  in  the  Teat,  in  Order  of  First  Appearance 
Page  Section  Symbol  Definition 

1  1  The  Laplacian  operator  with  reepect  to 

coordinates  (k); 

•  <51'/3ocw1  -t 

Atomic  number;  for  carbon,  6# 


1 

1 

2 

5 

4 
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2 


Atomic  orbital  for  electron  in  state  j 

with  coordinates  (k)  - K  or 

-  p(n^rk)S(tiT.i^(|j)o^(ei\  \)/\,  where 

P(x^^|r  )  is  the  radial  distribution  ampli¬ 
tude; 

is  the  spherical  harmonic  part; 
iB  the  spin  factor; 

and  the  quantum  numbers  have  the  usual 
significance* 
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Appendix  II 


W(1)=W (rx) 

V'>  ,  , 

T^nl+rilf  r)  =  -p^jdi 
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The  energy  parameter  of  the  equation  for 
the  J.*'*1  orb  ita  1  • 

The  Slater  exchange • 

r  P(n /  \r  )/r  •  _ 

dr*P(H/.  \  r^n'/jr)^  ♦  r^^- 


Lldr.P(n^|] 

J  *P(nT|r ' 
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r  *P(nT|r)/r^ 


Central  difference  op<££or • 
Associated  Legendre  polyncaials* 
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1*  Introduction-— A  number  of  approximation  methods  havd  been  applied  to  the 
many— body  problem  in  quantum  mechanics*  The  moat  accurate  one— electron 
method  is  that  of  Hartree  and  Pock  (v*  Appendix  I),  which  includes  the 
Pauli  principle  in  its  framework,  and  as  a  mathematical  consequence 
supposes  the  individual  electrons  to  move  in  potential  fields  which  are 
not  the  same  for  all*  Work  in  the  band  theory  of  solids  and  problems 
in  molecular  binding,  etc*,  requires  wave  functions  arising  from  a  poten¬ 
tial  field  which  is  the  same  for  all  electrons— a  "unipotential"— if 
calculations  are  to  be  practicable*  J  •  C*  Slater  has  proposed  ad  hoc  an 
averaging  process  to  be  applied  to  the  Hartree-Pock  field^  in  hopes  of 
retaining  the  accuracy  of  the  Hartree-Fock  method  while  accruing  unipoten¬ 
tial  simplicity* 

In  this  thesis,  the  Hartree-Fock  field  is  evaluated  for  the  configu- 

2  5  5  o 

ration  Is  2s2p^t  S  of  carbon  from  wave  functions  given  by  A*  Jucys,  and 

is  tabulated  in  such  a  manner  as  to  facilitate  calculation  of  unipotentials. 

Slater's  and  others*  The  first  phase  of  a  self-consistent  calculation 

is  then  carried  out  using  Slater's  unipotential;  the  set  of  wave  functions 

for  the  valence  electrons  which  obtains  differs  little  from  Jucys* • 

The  Hartree-Fock  equations  have  been  modified  for  tractability  by 
4  1 

Dirac  and  Slater;  for  an  atomic  system  the  modified  equations  reads 


1.1  u  >  ]  K/0=  ;M) 

(cf*  Appendix  I),  where  (1)  represents  the  codrdinates,  including  spin, 
of  the  electron  in  the  JLth  state,  and  (2)  the  corresponding  umbral  variables 
of  integration;  the  integral  sign  implies  summation  over  spin*  In  the 
literature,  the  first  summation  is  called  the  Coulomb  term  of  the  poten¬ 
tial,  while  the  second  is  usually  designated  the  exchange  term*  Hartree 


1.  J.  0*  Slater,  Phys*  Rev*,  81,  585,  '51. 

2*  A*  Jucys,  J.  Phys*  USSR,  XI,  #1,  '47* 

5*  R*  T*  Sharp  and  G*  K*  Horton,  Fhys*  Rev*,  90,  5^7,  '5% 
4*  P*  A*  M.  Dirac,  Proc «  Camb*  Phil*  Soc *,  26,  *50. 
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unite  are  used  here  and  throughout  this  thesis. 

Only  the  exchange  terra  of  1.1  is  not  the  same  for  all  electrons* 

In  the  case  of  the  Is  electrons  of  quadrivalent  carbon,  this  non-uniform¬ 
ity  is  not  an  objection  in  the  sense  of  paragraph  one  above,  for  the  K 
electrons  are  not  involved  in  interaction  calculations*  (Moreover,  self- 
consistent  calculations  for  the  beryllium  atom  and  ions  and  for  the  carbon 
atom  have  shown  the  core  electrons  to  be  insensitive  to  the  presence  of 

6  7 

other  electrons  and  to  exchange  effects;  so  that  for  this  work,  the  radial 

p 

part  of  the  Is  wave  functions  may  be  taken  as  that  given  by  either  Jucys 
£ 

or  Torrancel  the  difference  is  at  most  *002*)  For  the  valence  electrons, 

A 

however,  the  exchange  must  be  rendered  uniform;  this  was  accomplished 
by  Slater's  method*  While  the  above  procedure  amounts  to  calculating  the 
Is  and  2s  radial  parts  froa  different  potentials,  the  objection  that  they 
will  no  longer  be  orthogonal  is  not  serious  in  practice;  the  error  intro¬ 
duced  by  ignoring  the  defect  from  orthogonality  (as  it  turns  out,  a  small 
defect)  is  probably  obscured  by  the  intrinsic  inaccuracy  of  the  approximate 
theory  under  consideration!  moreover,  one  of  the  Is  orbitals  is  orthogonal 
to  the  2s  orbital  because  of  opposed  spin* 

The  exchange  appropriate  to  the  valence  electrons  is  shown  in  Appen¬ 
dix  I  to  bet 

>  w(D  -  r  i:  (*•>  k  'u-(x) '“-»«)  / 21  M<'i) 

where  the  1  summation  is  over  valence  electrons  only* 

In  order  that  equation  1*1  and  its  counterpart  with  exchange  1*2 
shall  be  separable,  the  potential  must  be  assumed  spherically  symmetrical* 

5*  D.Hartree,  Proc *  Roy*  Soc.,  14 1,  282, 

6*  D*  R*  and  W*  Hartree,  Proc  *  Roy*  Soc*,  A 15  A,  ’5^* 

7*C-Torrance,  Phys *  Rev*,  46,  5®8,  ’54* 
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Fortuitously,  this  condition  is  met  by  the  aggregate  integrands  of  the 


Coulomb  term  for  quadrivalent  carbon;  where  angular  dependence  arose  in 

ecA 

the  exchange  it  was  suppressed  by  an  averaging  process  effect^ by  multiply¬ 
ing  by  d^d^sin#,  and  integrating  over  all  angles* 

The  assumption  «f  a  spherically  symmetrical  potential  gives  rise  to 
the\»ual  spherical  harmonics  for  the  angular  parte  of  the  u^(>);  only 
the  radial  equation  need  be  solved.  This  reads: 

,  \  -  a.  [E^  -  v] 

where  V(r)  includes  the  appropriate  exchange  term*  Under  the  transform¬ 
ation  P(n£\r)  *  rRjj (iT)t  the  radial  equation  becomes 

1.5  ^  -  x  Ce^  -  v)]  pm  k)  . 

2o  Evaluation  of  the  Coulomb  and  Exchange  Terms— To  facilitate  manipulation 
of  the  integrals  appearing  in  the  Coulomb  and  exchange  terms,  indices 
were  assigned  to  the  states  occupied  by  the  electrons  in  the  conf iguration 


p  55 

l»<£2s2p':  3  as  follows. 

Index 

1 

2 

4 

5 

4 

5 

6 


Electronic  State  nlnnjm 
1  0  0  1 

1  0  0-3 

2  0  0 
2  1  1 
2  10 
2  1-1 


The  integrals 


2.1 


and 


2.2 


■  -  j  A  ^  ^  (*-)  > 


V*- 

called  respectively  the  Coulomb  and  exchange  integrals,  were  reduced  to 
single  integrations  over  the  variable  r^.  The  details  are  giv^n  in 
Appendix  II,  where  1,^  is  worked  out,  since  it  exemplifies  all  the 

*  kv 

pecularities  encountered,  including  angular  dependence. 
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Expressions  follow  for  the  various  integrals,  averaging  over  angles 


hav^ing  been  effected  where  necessary*  Integrals  of  the  form  X.*  I  ’ 
j  ^  2,  vanish  because  of  opposed  spin* 

The  Coulomb  term,  as  evaluated  by  the  method  of  Appendix  II  is 

5  I(r)  -fl,  -  ^  r2r^(l8,l«;r)  +  Y0(28,2B;r)  +  5Y0(2p,2pir)| 

The  exchange  integrals  are  given  by: 


Px(2s\r)Y0(2s,2s;r)/47ir‘i  • 

P"(2p\r)(Y0(2p,2p;r)  «•  ^Yj2p,2p;r)l  /hr,T) 

“  ^ ; 

F^(2pW)(jo(2p,2pJr)  A.  ^.Yj2p,2psr)7/A^r3,< 
I=  P(ls\r)P( 2»\r)Y  (le,2»;r)/4nr'i  ■ 

1  ^  O  ) 

I  -  P(lslr)P(2p'r)r,(lB,2p;r)/12nr 


=  I  ■  I  • 

15  <t  ■> 


1^*  P(2s|r)«2p|r)£Y  ( 28, 2p; r )J/l2ir  r  3 

=  L?  . 

5Pi(2p'r)Yi.(2p,2p;r)/100<r  r  3; 

”  £*»  =  I,  • 

ft 

The  foregoing  were  evaluated  numerically  from  Jucys'  wave  func- 

2 

tions  by  the  following  method*  Let  f(r)  be  an  integrand  tabulated  at 
intervals  h  of  the  argument;  then  integrals  of  the  form 


I  f(r)dr  were 

l  ,  44 

evaluated  by  iimpson's  rule,  where  k  is  an  integer.  The  integral  I  f(r) 

/  4 

was  then  calculated  by  Cotes'  three-eighths  rule,  whereupon  {  f(r)dr  = 


dr 


( 

f(r)dr  -  f 

Jo 


f(r)dr  was  found*  To  complete  evaluation  of  the  integral 

up  to  all  tabular  values  of  r,  integrals  of  the  form  lf(r)dr  were 

conputed  by  Simpson's  rule*  - 

<jre 

I(r)  as  given  by  2*5  is  tabulated  in  Table  I,  as.the  combinations 
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which  land  then  selves  to  calculation  of  the  unipotential  in  this  and 

5 

other  problems* 

As  shown  in  Appendix  I, 

2.5  W(r)  -  (jS(r)  4  0(r)]/^§1(2.|r)  4  ?p\2p  lr)J  I 

this  expression  is  also  tablulated  in  Table  I* 

Plots  of  I(r)  and  W(r)  vs*  r  are  given  in  Fig*  I* 

5*  Solutions  of  the  Radial  Equation  for  Small  and  Large  r — Equation  1*5  may 
be  written 

%1  d2P(n£\r)/dr*'=  \l  (1*1)  -28/r  4  2(1  -  W)  4  2^^  P(n^\r)  . 

(<jl)  In  order  to  obtain  a  working  approximation  to  the  solutions  of  5*1 
in  the  vicinity  of  the  origin,  a  cubic  curve  was  fitted  to  four  consee- 
utive  values  of  ( I  -  W),  starting  with  the  second  from  the  origin,  by 


means  of  Lagrange's  interpolation  formula*  For  F(r)  -  (I-W)  =  a  +br 
**  ar  4  dr,  Lagrange's  interpolation  formala  yields: 


c  =  — V  \  ]  ~  FJr,+  K*  -  *[Fi<*\*  V  *  p,(r,+ITIV^ 

d=  -£(V  f,  )  -  i(F^-  rz), 

where  F(fr^)  s  Fi#  and  r  =  *02,  r^*  *04,  etc* 

Equation  nay  now  be  written  in  the  form 

5*2  d'tydr' 4  (<*/r"4-  £/r  +1/  +  5r  4  trlv  fr3)P  *  0;  (r  _<  0.08  ) 

here  *s-/(/4l),  [3  =  28  =  12,  and  y  is  a  augmented  by  2^^^  •  This 

equation  was  solved  by  Frobeniue'  method;  writing 

8 *H .Jeffreys  and  Jeffreys,  Methods  of  Mathematical  Physics,  p«  527* 

1  Camb^d^e  Univ'e^sit^  fVe^s; 
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p  =  21  a^r^f  ( 0,  since  P(n^|o)  =  0), 
and  substituting  into  5*2  yields  recursion  formulae  for  the  a„by  means 
of  which  two  ordinates  for  purposes  of  starting  numerical  solution*  out¬ 
ward  from  the  origin  were  calculated*  In  Table  II  are  given  the  a^  from 
which  starting  ordinates  for  the  final  solutions  were  calculated;  it  is 
indicative  of  the  inse^nsitivity  of  the  solution  to  changes  in 
the  vicinity  of  the  origin)  that  while  the  aychanged  with  changes  in\£x^| 
the  ordinates  P(n^*02)  and  P(n-£\,o4)  did  not* 

•ince  the  series  5*5  is  certainly  convergent  as  a  consequence  of 
the  regularity  of  the  differential  equation  of  which  it  is  a  solution^, 
and,  as  it  turns  out,  alternating,  truncation  error  is  smaller  in 
absolute  value  than  the  first  neglected  or  last  counted  term* 

(b)  By  Table  I,  for  r2^*0,  l(r)  =  6/r,  and  W(r)  *  1/r,  so  that 
5*1  may  be  written 


d  P(n£\r)/drL-  )/r-2/r  2  Ie^  P(n£|r),  r}_6* 0; 

this  equation  is  soluble  in  terms  of  confluent  hypergeoraetric  functions: 

10 

asymptotic  expansions  valid  for  large  r  are 

P(2.|r)  ^  A.-r/Ir*  [  \  +  "  )  •  1  =  > 

P(2pW)-A  r*f  1*  \  .  K  - 

<?a.r/K  J  /a. 

Table  II  gives  numerical  values  for  the  various  constants*  These  express¬ 
ions  were  used  to  calculate  starting  ordinates  for  the  inward  numerical 
solutions  dealt  with  in  section  4  below* 

4.  Numerical  Solution  of  the  Radial  Equation — The  physically  significant 
solutions  of  5*1  result  from  choosing  a  value  of  the  energy  eigenvalue 
'l»n{  which  yields  a  Veil-behaved"  solution,  i*e.,  a  solution  which  approach¬ 
es  zero  with  a  slope  decreasing  in  absolute  value  for  large  r*  The 


9*  E*  Whittaker  and  0*  Watson,  Moddrn  Analysis,  p*199» 
10*  Whittaker  and  Watson,  Ibidem,  p*  557  ff • 

(  O*-' '  '/fiiTS  t  ly  SS  t  <S,  134-4) 
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asymptotic  expansions  certainly  meet  this  last  requirement*  Eval¬ 
uation  of  the  eigenvalue  was  accomplished  by  the  method  used  by  the 
Hartrees:^  for  a  given  estimate  of  numerical  solutions  of  5*1  were 

computed  from  the  origin  outward,  and  from  r  =  6*5  inward,  using  for  a 
start  ordinates  computed  by  the  methods  of  section  if  the  solutions 
could  be  made  to  merge  smoothly  (in  the  vicinity  of  r  :  5*5)>  then  the 
eigenvalue  was  adjudged  correct*  "Smoothness "  was  judged  on  this  cri¬ 
terion:  ordinates  were  calculated  from  each  solution  for  r  =  5*5  and 
r  s  4*0;  if  the  respective  ordinates  could  be  rendered  equal  (within 
estimated  error)  for  both  of  these  abscissae  by  the  same  numerical 
factor,  then  junction  was  considered  smooth*  Tables  IV  and  V  show  that 
confluence  within  a  few  units  in  the  fifth  decimal  place  was  obtained*. 

(The  choice  of  r  5  5*5  as  the  abscissa  of  confluence  was  dictated 
by  the  intervals  in  the  tables  of  I(r)  and  W(r)j  fortuitously,  it  deve¬ 
loped  that  the  estimated  cumulative  errors  due  to  tnmcaticwi  of  a  differ¬ 
ence  series  used  in  the  solution  were  of  the  same  sign  and  order  of 
magnitude  there • ) 

To  implement  numerical  solution,  equation  was  written 
1  d^Ptn^rJ/dr  =  gnJ?(r)P(n/\r), 

with  gn^(r)  =  l{tsl)/r~-  2&/r  +  2(1  -  W)  2^/].  The  function  g^was 
calculated  using  the  functions  appearing  in  Table  I,  and  tabulated;  the 
final  values  of  gn^  are  tabjuilated  in  Tables  IV  and  V,  incidentally  to 
dee*onst rating  the  method  of  solution  of  4.1  for  the  2s  and  2p  cases 
respectively* 

The  method  of  solving  4.1  used  is  described  by  Jeffreys  and  Jeff- 
12 

re7**  By  definition  of  the  second  central  difference 

11.  Hartree  and  Hartree,  Proc*  Roy*  Soc*,  Al^4,  588,  *^6 

12*  Jeffreys  and  Jeffreys,  oj>.  cit*,  p*  267  ff * 
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o^P(r)  =  P(rfh)  -  2P(r)  4  P(r-h) 

=  (eM  -  2  4  ©“hD)P(r) 

-  £(co%h  hD  -  l)P(r) 

*  |D  p^), 

so  that,  symbolically, 

4*2  hD/2  »  ■s/'»A_,(-^r) 

Further, 

S"P  =  hVp  +  b*1"'’  *M>f  -  h*Dj  P 

«  h*gP  t  C(  V'"*i  %hn/ihD)2  -  lj  h^gP 
by  equation  4.1.  Using  4*2,  one  obtains 

4.5  Vp  =  tf“gP  4  h''||S/5,o/,-,(iS  . 

The  operator  in  brackets  may  be  expanded  in  a  power  series;  when  this 
is  done,  equation  4*5  reads 

4.4  5Vh'  -  gP  -*•  ^(gP)/^  -  S4(gP)/240-  ...  # 

(Equation  4.4  is  derived  in  a  more  straightf orward  but  laborious  manner 

by  £-.  Milne}?  .) 

Consider  now  a  function  denoted  by  V’tgP)  whose  second  differences 
are  gP;  in  terms  of  such  a  function  4.4  can  be  written 

4.5  SVh1'-  ^[^gP)  ♦  gP/12  -  S"(gP)/240  ♦  •••], 

fron  which  is  apparent  that  the  functions 

4.6  P/h*"  and  [_5"^(gP)  4  gp/i2  -  S*(gP)/240  4  »**~\ 

have  the  same  second  differences.  The  above  definition  of  5  -l(gp)  does 
not  determine  it  uniquely,  for  any  function  of  the  form  A4Br,  whose 

15.  £.  Milne,  Numerical  Solutions  of  Differential  Equations,  pp.85#  84. 

(wMevj  is>b3.^) 
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second  differences  vanish,  may  be  added  to  a  given  S'KgP)  an<*  the  second 
differences  of  the  sum  will  still  be  gP.  In  particular,  A  and  B  can  be 
so  chosen  as  to  render  the  functions  4.6  equal  for  two  adjacent  tabular 
values  of  P;  then  those  functions  will  be  equal  for  all  tabular  values 
of  the  argument,  wherefore 

4.7  P/lT -  ^(gP)  *  gP/12  -^*(gP)/240  -V  • 

In  practice,  4,7  may  be  truncated  at  the  second  term* 

4.7a  P/h^  *"*■( gp)  +  gP/12  . 

The  «*ror  accumulated  over  a  range  of  integration  due  to  truncation 
can  be  estimated  by  summing  hL5'’(gP)/240  over  that  rangej  this  was  effect¬ 
ed  by  differencing  the  leading  and  final  differences  of  gP  in  that  range 
(with  due  regard  for  sign)  and  multiplying  by  hV240. 

In  Table  III,  the  equation  d^P/dr* =  P  is  soVved  and  compared  with 
the  ascending  solution  P  =  e  to  demonstrate  the  accuracy  of  the  method. 

Equation  4.7a  was  applied  to  the  soluticxi  of  4.1  as  follows.  Values 
of  P  for  r=.02  and  .04  were  calculated  by  the  method  of  sec.  5  and  entered 
in  the  P  column  (v.  Table  IV t  the  starting  values  are  enclosed  in  paren¬ 
theses).  (The  knowA value  P(n^Jo)  =  0  could  not  be  used  because  £^(0)  i® 
indeterminate.)  Fran  the  entries  under  P,  the  corresponding  values  of 
P/h1-,  gP,  and  by  means  of  4.7a,  5  *{gP)  were  computed.  The  dirierence 
of  the  last-mentioned  was  then  entered  in  the  S~'(gP)  column,  and  added 
to  (gP)#o4  yi®ld  the  next  value  of  S~'(gP),  which  was  in  turn  added 
to  V'tgP)  to  give  V'TgP)  ^  • 

To  proceed  with  the  solutiwi,  the  value  of  (gP)  was  estimated 

and  (P/h’")  calculated  therefrom  by  4.7a;  this  was  then  used  to  cotn- 

Oi, 

put©  a  value  of  (gP)  which  was  compared  with  the  estimated  value. 

■Of, 

If  different,  the  new  value  was  used  to  calculate  another  value  of  P/h^ 
and  this  process  was  repeated  until  no  further  difference  obtained. 
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The  new  value  of  gP  was  then  added  to  the  nearest  S~^gP)  to  yield  the 
next  value  of  5  '(gP),  which  was  then  added  to  the  nearest  5  l(gP)  to  get 
the  next  entry  under  $  (gP),  etc*  The  inward  solution  was  conducted 
in  a  exactly  similar  manner* 

At  a  change  of  he  tabular  interval,  it  was  often  necessary  to 
interpolate  a  value  of  P  halfway  between  two  know  values  (e*g.,  P(2^>25) 
Table  IV);  this  was  done  by  calculating  5  (gP)  on  the  new  interval 
for  the  known  values,  then  calculating  the  required  ordinate  from 

P(r)  =  (S'lV)r-h+S'VU|  /  C^+iVr'l 
fo  demonstrate  this,  the  numerator  of  the  right  hand  side  may  be  written, 
with  the  aid  of  4*7,  as 

S'V(*P)  ■*  =  (p(rt-h)  4  Ptr-h^/h1"-  CgP)  +  (gP)  1/12 

•  S"p/hXf  2P/h>*-  %\gP)/12  -  gP/i6,  by  4.4 

*  (2/h%  5g/6)P,  Q.E.D* 

Within  Tables  IV  and  V  are  shown  also  the  details  of  merging  inward 
and  outward  solutions* 

Table  VI  gives  normalized  values  of  the  P(n^lr)  and  their  eigenval¬ 
ues.  Table  VII  compares  the  eigenvalues  obtained  with  those  found  by 

2  14 

Jueys  and  Torrance* 

Pig.  II  is  a  plot  of  the  P(n£\r)  vs.  r  compared  with  the  correspond¬ 
ing  functions  obtained  by  Jueys* 

E. 

H  Levy  and  Baggot  point  out  that  to  yield  a  solution  correct  to  three 

A 

significant  figures,  say,  the  numerical  work  must  be  carried  to  five. 

All  extra  figure  was  carried  through^  in  the  calculation  of  o  (gP)  and 
gP  to  give  rounding-off  errors  a  place  to  accumulate  harmlessly^ 

14.  C.  Tor  ranee,  16c*  primo  cit.,  p*  5^9  • 

15*  H.  Levy  and  E*  A*  Baggott,  Numerical  Studies  in  Differential 

Equations,  vol.  I,  p.112.  (o ove<  V\»YdV>cx*T>ov\s  Nev  YcnrV;,  *35o) 

16.  H.  Jeffreys  and  B*  Jeffreys,  op*  cit*,  p277j  also  B.  Hartree, 

Numerical  Analysis,  pp.  85-84.  I  C Navev^ov*  Oxrfovd, 

\S52.  ) 
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TABLE  I 


r 

I(r) 

B(r) 

C(r) 

W(r) 

0.00 

14.77 

_ — - 

— 

0.02 

14.75 

27.38 

0.1414 

0.9428 

0.04 

14.44 

89.43 

0.5760 

3.922 

0.06 

14.13 

70.76 

1.002 

4.007 

0.08 

13.71 

53.08 

1.520 

3.886 

0.10 

13.30 

38;61 

1.974 

5.646 

0.12 

12.84 

28.66 

2  .  c68 

3.523 

0.14 

12.40 

20.00 

2.659 

3.234 

0.16 

11.95 

14.59 

2.872 

3.098 

0,18 

11.48 

9.839 

5.060 

2.795 

0.20 

11.10 

6.999 

3.167 

2.637 

0.25 

10.15 

2.524 

3.224 

2.096 

0.-3O 

9.355 

1.531 

3.230 

2.014 

0.35 

8.664 

+0.08489 

5.116 

1.383 

0.40 

8.088 

-0.00394 

3.028 

1.243 

0.45 

7.592 

+0.1065 

2.924 

1.170 

0,50 

7.176 

0.2829 

2.828 

1.140 

0.55 

6.810 

0.4478 

2.730 

1.120 

0.60 

6.498 

0.5945 

2,641 

l.H3gQ 

0.65 

6.218 

0.7175 

2-529 

1.114* 

0.70 

5.974 

0.7765 

2.419 

1.106" 

0.75 

5.749 

0.8145 

2.294 

1.101 

0.80 

5.550 

0.8324 

2.163 

1.098 

0.90 

5.178 

0.8027 

1.881 

1.081 

1.0 

4.839 

0.7257 

1.595 

1.057 

1.1 

4.528 

0.6596 

1.359 

1.059 

1.2 

4,245g 

0.5201 

1.080 

0.9840 

1.3 

3.997g 

0.4221 

0.8729 

0.9439 

1.4 

3.808q 

0.3372 

0.7073 

0.9115 

1.5 

3.608o 

0.2641 

0.5632 

0.8704 

1.6 

3.4504. 

0.2056 

0.4527 

0.8376 

1.7 

3.277., 

0.1579 

0.3550 

0.7954 

1.8 

3.143*; 

0.1211 

0.2821 

0.7622 

1.9 

2.992® 

0.09168 

0.2205 

0.7225 

2.0 

2.8795 

0.06966 

0.1740 

0.6922 

2.1 

2.754a 

0.05252 

0.1353 

0.6555 

2.2 

2.4>500 

0.03953 

0.1068 

0.6285 

2.  3 

2.5316 

0.02962 

0.08290 

0.5945 

2.4 

2.4501 

0.02227- 

0.06538 

0.5710 

2.6 

2.2754 

0.01252 

0.04017 

0.5305 

2.8 

2.1219 

0.00687 

0.02455 

0.4710 

3.0 

1.9864 

0.00395 

0.01517 

0.4377 

5.2 

1.8660 

0.00221 

0.009535 

0.4034 

3.4 

1.7590 

0.00106 

0.00565^ 

0.3530 

3.6 

1.6620 

0.000697 

0.003626 

0.3455 

3.8 

I.5765 

0.000595 

0.00227e 

0.5  18 

4.0 

I.4983 

0.00022 

0.00143. 

0.3005 

4.5 

1.5325 

0. 00005® 

O.OOO453 

0.2576 

5.0 

1.1996 

O.OOOOI4 

0.00014g 

0.2272 

5.5 

1.090q 

O.OOOOO4 

0.00004 

0.2041 

6.0 

l.OOOo 

0.00000 

O.OOOOI5 

0.1667 

subscript  numb ora  are  only  partially  significant. 
W(r)  -  1/r,  I(r)  m  6/r,  for  r>  6.0. 
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Figure  i 

COULOMB  (I)  AND  EXCHANGE  (W  )  TERMS 

(HARTREE  UNITS) 
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table  ii 


P(nllr)  s  ArC(a0  t  a,r  *  a^r2  f  ***)  9  .04 


nl  A  c 

2b  -6.08804  1 

2p  5.O25I5  2 


O  6  9 

1  -6  1.9896  -149.852  158O.595  _ 

1  -1/5  5.051  -75.0990  815.5515  -5171.58 


P(nl|r)  =  Ar^e  cr(l  -*>  a^/r  ag/r^  ■»  a^/r^  t  •••),  r  ^  6.0 
nl  A  b  c  a ,  aa 

2e  -5.528O  o79606ll  I.256I85  .0646194  -.00515757 

2p  1 .62085  (-)  1.047824  .9545584  1. 02157}  *547599'  .0154968 


table  III 

It 


Numerical  solution  of  P  =  P,  compared  with  ascending  solution  P  -  er. 


r 

S" 

8-'p 

P 

loop 

r 

<» 

error 

.00 

• 

99.916 

10.511 

(l.OOOO  ) 

(100.00) 

1.0000 

.0000 

.1 

110.428 

11.616 

(I.IO52  ) 

(110.52) 

1.1052 

0 

.2 

122.044 

12.858 

1.222 

122.15 

1.2214 

1 

•5 

154.882 

14.188 

1.550 

154.99 

1.5499 

0 

.4 

1  49.070 

15 .680 

1.492 

149.19 

1.4918 

1 

•5 

164.750 

17.529 

1.6  4o 

164.89 

1.6487 

2 

.6 

182.079 

19.151 

1.822 

182.25 

1.3225 

0 

•7 

201.250 

21.165 

2.014 

201.40 

2.158 

2 

.8 

222.595 

25.591 

2.226 

222.58 

2.2255 

5 

.9 

245.7 86 

25.851 

2.460 

245.99 

2.4596 

5 

1.0 

* 

Parenthesized 

271.657 

entries  are  "starting 

2.719 

;  entries'] 

271.86 

2.7185 

.0005 

«•  +  * 
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TABLE  IV 


Solution 

of  d*P/dr  = 

gP  for  2s  case,  unnormalized.  2S^ 

-1.5780. 

r 

&-(sp) 

S"(sp) 

(gP) 

25OOP 

g 

P 

0.02 

(275.157) 

(207.454) 

-(61.647) 

(270.00) 

-570.8020 

(CU0800) 

.04 

(482.571) 

154.582 

(55.052) 

(478.15) 

277.5820 

♦ 

(0.19126) 

.06 

656.955 

109.255 

45.127 

655.19 

1784720 

.25528 

.08 

746.208 

70.986 

58.269 

745.02 

128.7620 

.29721 

.10 

817.19* 

58.692 

52.294 

814.50 

99  .1220 

.5258O 

.12 

855.886 

♦  11.447 

27.245 

855.62 

79.7920 

.54145 

.14 

8^7.555 

-11.552 

22.779 

865.45 

65.8020 

.54617 

.16 

856.OOI 

50.572 

19.040 

854.41 

55.7120 

.54176 

.18 

825.629 

46.104 

15.752 

824.52 

47.7120 

•52975 

.20 

779.525 

59.021 

12.917 

778.449 

41.4920 

.>1158 

.22 

.24 

.26 

.24 

720.504 

65I.O7I 

575.505 

2602.20 

69.455 

-77.766 

-151.61 

10.412 

8.555 

-6.655 

-8.555 

719.64 

650.577 

572.752 

4 

10  P 

2601.51 

56.1720 

52.0520 

-28.9520 

-52.0520 

• 

•25 

.26 

2450.59 

2291.56 

-*•159.04 

7.45 

-6.65 

2449.97 

2291.01 

50.5120 

-28.9520 

.24500 

.20 

125.6282 

-27.0105 

-12.9166 

124.5518 

-41.4820 

.51158 

•25 

98.6177 

54.4568 

7.4265 

97.9988 

50.5120 

.24500 

.50 

64.1809 

58.2275 

<•5.7907 

65.865 

25*7420 

.15966 

•55 

♦25.9554 

59.4008 

-1.1755 

♦25.856 

18.1520 

♦  .06464 

.4o 

-15.4474 

58.9071 

♦0.4957 

-15.406 

14.7520 

-0.0555 15 

.45 

52.5545 

57.507 6 

1.5995 

52.221 

12.2520 

.150550 

.50 

89.6621 

54.9921 

2.5155 

89.459 

IO.552O 

.22567 

•55 

124.6542 

-52.2556 

2.7565 

124.42 

8.8620 

.51105 

.60 

T  156.8898 

2.9965 

-156.64 

-7.6520 

•<59160 

Sign  is  indicated  only  at  the  head 
change. 

and  foot  of 

a  column,  and  at  a 

TABLE 

IV  Co  cat.) 

r 

5~v(gP) 

*“(gP) 

gP 

400P 

g 

P 

0.60 

-156.8898 

-29.2591 

2.9965 

-156.64 

-7.6520 

-.59160 

.65 

186.1289 

26.1588 

5.1005 

185.87 

6.6720 

.46468 

.70 

212.2677 

25.0488 

5.O9OO 

212.010 

5.85OO 

.55005 

•75 

255.5165 

5.0124 

255.07 

5.1260 

-.58768 

•80 

-20.0564 

-255.5529 

2.8815 

-255.115 

-4 .5 180 

loop 

.70 

-55.2605 

-IO.758I 

5.0901 

-55*0028 

-5.85OO 

-.550025 

.80 

64.0184 

7.87 66 

2.88 15 

65.7785 

4.5 180 

.657785 

•90 

71.8950 

5*5240 

2.5526 

71.682 

5.5610 

.71682 

1.00 

77.2190 

5.125* 

2.2009 

77.056 

2.8570 

•770>£ 

1.1 

8O.542I 

1.8586 

1.2645 

80.257 

1.5760 

.80257 

1.2 

82.2007 

-0.5009 

1.5577 

82.071 

1.8980 

.82071 

1.5 

82.5016 

♦0.9721 

1.2750 

82.596 

1.5450 

.82596 

1.4 

81.5295 

1.9478 

0.9757 

81.448 

1.1980 

.81448 

1.5 

79.5817 

2.7000 

.7522 

79.519 

9.9460 

•79519 

1.6 

76.8817 

5.2548 

•5548 

76.857 

.6960 

.76857 

1.7 

75.6469 

5.6154 

.5806 

75.615 

.5170 

.7*615 

1.8 

70.0515 

5.8456 

.2282 

70.012 

.5260 

.70012 

1.9 

66.1879 

5*9746 

.1510 

66.177 

.1980 

.66177 

2.0 

62.2155 

4.0058 

♦.0292 

62.211 

-.0470 

.62211 

2.1 

58.2095 

5*9685 

-.0555 

58.212 

♦  .0610 

.58212 

2.2 

54.2412 

5*8782 

.0901 

54.249 

.01660 

•54249 

2.5 

50.5650 

♦5*7598 

.1184 

50.575 

.2550 

.50575 

2.4 

-46.6052 

-.1566 

-46.616 

♦  .5560 

-.46616 

Sign  ie  indicated  only  at  the  head  and  foot  of  a  column,  and  at  a  change. 
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r  r 


r 

$"(gp) 

TABLE 

gP 

IV  (Cont.) 

25P 

E 

P 

2.2 

-15*5548 

-.0901 

-15.5625 

♦0.1660 

-0.54249 

2.4 

11.6410 

♦  1.9158 

•  I566 

11.6540 

.5560 

.46616 

2.6 

9.8858 

1.7572 

.1875 

9.899 

.4750 

.59596 

2.8 

8.5159 

1.5699 

.1  980 

8.550 

*594? 

.55520 

5.0 

6.9420 

1.5719 

.1879 

6.958 

•6750 

.27852 

5.2 

5*7580 

1.1840 

.1759 

5.772 

.7550 

.25088 

5.-4 

4.7479 

1.0101 

.1659 

4.762 

.8606 

.19048 

5.6 

5.9017 

0.8462 

.1577 

5.915 

.8797 

.15652 

5.8 

5.1952 

.7085 

.1195 

5.205 

.9511 

.12812 

4.0 

-2.6040 

.5892 

-.1018 

-2.612 

.9740 

-.10448 

5.4 

-19.054 

-.164 

loop 

-19.048 

t.8606 

-.19048 

5*5 

17.261 

1.775 

C\J 

in 

rH 

• 

17.274 

.8814 

.17274 

5.6 

-15.641 

1.021 

-.158 

-15.652 

.8797 

-.15652 

6.5 

-0.496244 

-0.16 18 15 

4p 

-0.509^28 

1.2698 

-0.127452 

6.0 

•875775 

-0.577551 

.279155 

.897056 

1.2447 

-  .224259 

5*5 

5.0 

4.5 

4.0 

5.5 

•656666 

1.55044 

I.II517 

2.64561 

1.87572 

4.52155 

5*1021(6 

7.62581 

-4.99752 

-12.62115 

Outward  Sol^.  " 

.45850 

.76055 

1.22676 

1.89484 

-2.8550 

Inward  Sol'ii. 

1.5686 

2.7090 

4.6256 

7.7817 

-12.857 

Conv's'n. 

1.1692 

I.I25O 

I.06I5 

0.9740 

.8814 

Factor  Conv.  I.  S. 

5.5 

-.17274 

-12.857 

.01545548 

17274 

4.0 

-.10448 

-7.7817 

“  • 

10455 

Normalizing  factor  for  this  solution*  0,990624 
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TABLE  V 


Solution  of  dP/dr=gP,  for  2p  case,  with  2E  -  -*9108;  unnormalized. 


r 

S'^gp) 

&-’(gP) 

gP 

2500P 

g 

P 

0.02 

4.262 

I5.OI5 

8.857 

5.000 

4428.5508 

0.00200 

.04 

19.275 

22.758 

7.745 

19.920 

971.9508 

.007968 

.06 

42.055 

29.172 

6.414 

42.567 

576.7164 

.017027 

.08 

71.205 

54.418 

5.246 

71.642 

I85.070 8 

.028657 

.10 

IO5.625 

58 .666 

4.248 

IO5.98 

100.2108 

.042592 

•  12 

144.289 

42.045 

5079 

144.57 

58.  2968 

.057828 

.14 

I86.554 

44.699 

2.654 

I86.56 

55.5716 

.074624 

•  16 

251.055 

46.710 

2.011 

251.20 

21.7458 

.091480 

.18 

277.745 

48.194 

1.484 

277.87 

15.5492 

.11148 

.20 

525.957 

49.218 

1.024 

526.022 

7.8508 

.150409 

.22 

575*155 

49.891 

0.675 

575.211 

4.4851 

.I5OO84 

•  24 

425.046 

50.255 

♦0.544 

425.075 

2.02502 

.170050 

.26 

475.281 

-0.006 

475*281 

10^p 

-0.0554 

.190112 

.24 

1700.27 

IOO.54 

♦0.54 

I700.50 

♦2.02502 

.17005 

•25 

1800.61 

IOO.52 

♦  .18 

1800.62 

♦1.0208 

•  180062 

.26 

1901.12 

-o.oo6 

1901.12 

400P 

-O.O554 

•190112 

.20 

52.0785 

19.9512 

I.O258 

52.I656 

7.850 8 

.150409 

•25 

72.0095 

20.H50 

♦O.I858 

72.0248 

♦  1.0208 

.180062 

•5° 

92.1245 

19.6115 

-O.5055 

92.085 

-2.1870 

.25021 

•55 

111.756 

18  16 

.696 

111.678 

2.4927 

.27920 

.40 

150.652 

17.970 

.946 

150.57 

2.8992 

.52645 

♦45 

148.622 

16.840 

1.150 

148.55 

5.0427 

.57155 

.50 

I65.462 

15.592 

1.248 

165.56 

5.0192 

.41540 

•55 

18 1.054 

14.272 

1.520 

I8O.94 

2.9176 

.45255 

.60 

195.526 

1.549 

195.21 

2.7656 

.48805 

Sign 

is  indicated  only  at 

the  head  and 

foot  of  a  column,  and  at 

changes. 
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TABLE 


r 

^"(gP) 

s"(gp) 

gP 

0.60 

195.526 

12.925 

-1.549 

.65 

208.249 

11.567 

1.556 

.70 

219.816 

10.240 

1.527 

•75 

250.056 

8.954 

1.286 

.80 

259.010 

-O.659 

(Cont.) 


400P 

g 

P 

195.21 

-2.7656 

0.48805 

208.14 

2.6O55 

.52055 

219.705 

2.4156 

•54926 

229.95 

2.2576 

.57488 

258.955 

-1.1052 

.59759 

.70 

55.0569 

.80 

59.7957 

•90 

65.8915 

1.0 

66.8670 

1.1 

68.8246 

1.2 

70.5761 

1*5 

71.1006 

1.4 

71.0942 

1.5 

70.4876 

1.6 

69.5708 

1.7 

67.8505 

1.8 

65.9964 

1.9 

65.8945 

2.0 

61.5954 

2.1 

59.1606 

2.2 

56.6575 

2.^5 

54.0646 

2.4 

51.4625 

4.7568 

4.0978 

2.9755 

1.9576 

1.5515 

♦0.7245 
-0.006 A 
.6066 
1.1168 
I.5205 
1.8541 
2.1021 
2.5OO9 
2.4528 

2.5251 

-2„5729 

-2.6021 


-1.5268 

O.659O 

1.1225 

1.0179 

0.4061 

.8270 

.7509 

.6002 

•5102 

.4055 

.2480 
.1988 
.1519 
.0905 
.0498 
-.  0292 
♦.0082 


100P5 

54.9265 

59.7588 

65.798 

66.782 

68.791 

70.507 

71.040 

71.044 

70.445 

69.557 

67.825 

65.976 

65.878 

61.582 

59.155 

56.6554 

54.062 

51.4652 


-2.4156 
I.IO52 
1.759* 
1.5242 
°.5905 
1.1765 
1.0288 
0.8448 
.7245 
•5819 
.4922 
•5759 
•5112 
.2142 
.1527 
.0880 
-.0541 
♦  .0160 


0.549265 
.597588 
.65798 
.66782 
.68791 
.70507 
.71040 
.71044 
.70445 
.69557 
.67825 
.65976 
.65878 
.6 1582 
.59155 
.566554 

.54062 

•514652 
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TABLE  V  (Cant.) 


r 

S'XsP) 

S-’(gP) 

gP 

2.2 

14.1625 

-1.2974 

-.0498 

2.* 

12.8651 

1.2892 

♦.0082 

2.& 

11.5759 

1.2421 

•0471 

2.8 

10.5558 

1.1667 

.0754 

5.0 

9.1671 

1.0825 

.0844 

5.2 

8.0848 

0.9915 

.0910 

3.4 

7.0955 

.8872 

.1041 

3.6 

6.2065 

.7960 

.0912 

5*8 

5.4105 

-.7088 

.0872 

4.0 

4.7015 

.O8I5 

3.* 

28.400 

-1.8260 

.104 

5.5 

26.574 

-1.726 

.100 

5.6 

24.848 

.091 

6.5 

0.66579 

O.555I5 

0.10955 

6.0 

0.99694 

.49502 

.15987 

5.5 

1.48996 

I70717 

.21415 

5.o 

2. 19715 

1.00479 

.29762 

4.5 

5.20192 

1.40544 

.59865 

4.0 

4.60556 

1.90511 

.50167 

5.5 

6.51047 

.61950 

Outward  Sol'n.  Inward  Sol'n. 

3.5 

O.26582 

6.562I 

4.0 

.I8855 

4.6472 

25P 

g 

P 

14.1584 

-0.0880 

0.566554 

12.8658 

♦0.0160 

.514652 

II.58O 

.1017 

.46520 

10.540 

.1822 

.41560 

9.174 

.2500 

.56696 

8.092 

.2811 

.52568 

7.1022 

•5664 

.28409 

6.2159 

.5668 

.24856 

5.418 

.4024 

.21672 

4.7085 

.4518 

.I8855 

loop 

28.409 

.5664 

.28409 

26.582 

•5775 

.26582 

24.856 

.5668 

.24856 

4p 

0.672900 

0.6499 

0.168225 

1.01026 

.6550 

.252565 

1.5078 

.5681 

2.2219 

.5558 

5.2551 

.4929 

4.6472 

.4518 

6.5621 

.5775 

Conv’s'n  Factor  Converted  I.  S. 

.0405084  O.26582 

.18825 


Normalization  factor  for  thie  solution:  I.OOO5I5 
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TABLE  yl 


Normalized  solutions  of  the  2b  and  2$  radial  equations* 


r 

P(2s[r) 

P(2p(r) 

r 

P(2e|r) 

P(2pfr) 

0.00 

0.000 

0.000 

1.7 

-0.729 

0.678 

.02 

.107 

.002 

1.8 

.694 

.660 

*04 

.190 

.008 

1.9 

.656 

.659 

.06 

.251 

.017 

2.0 

.612 

.616 

.08 

•294 

♦029 

2.1 

•577 

.592 

.10 

;525 

.042 

2.2 

•557 

.587 

.12 

•558 

.058 

2.5 

•4 99 

.5^1 

.14 

.075 

2.4 

.462 

.515 

.16 

•559 

.095 

2.6 

•592 

.465 

.18 

•527 

•  111 

2.8 

•550 

.414 

0 

CM 

. 

.509 

.150 

5.0 

.27  6 

.587 

•25 

.245 

.180 

5.2 

.229 

•524 

.50 

•158 

.250 

5s* 

.189 

.284 

•55 

♦0.064 

.279 

5.6 

.155 

.249 

•4o 

-O.O55 

•527 

5.8 

.127 

.217 

.45 

.129 

•571 

4.0 

.104 

.188 

.50 

.222 

.414 

4.5 

.062 

•  151 

•55 

•508 

•  ^55 

5.0 

.056 

.090 

.60 

.588 

.488 

5.5 

.021 

.061 

.65 

.460 

•521 

6.0 

.012 

.041 

.70 

•525 

•550 

6.5 

.007 

.027 

.75 

•582 

•575 

7.0 

.004 

.018 

.80 

.652 

•59  8 

7.5 

.002 

.012 

.90 

.710 

.658 

8.0 

.001 

.008 

1.0 

.7  65 

.668 

9.0 

•OOO4 

.005 

1.1 

.795 

.688 

10.0 

.ooo1 

.001 

1.2 

.815 

.705 

11.0 

•ooo0 

.OOOi 

1*5 

.816 

.711 

12.0 

.0000 

.0000 

1.4 

.807 

.711 

1*5 

.788 

.705 

Asymptotic  expansions 

are  given  in 

1.6 

-0.761 

0.694 

Table 

II. 

table  vii 

Comparison  of  eigenvalues  arising  from  various  exchange  terms* 


Exchange  -21  “2Ki» 

None7  0.6425  1*2454 

Hantree-Fock2  O.9569  I.885I 

Slater  O.9JL08  1*5780 
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RADIAL  DISTRIBUTION  AMPLITUDES 
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5*  0 one lua lone— The  wave  functions  tabulated  in  Table  71  are  not  self- 


consistent  on  the  usual  criterion  that  the  integrals  jdrPv(n/\r)  eal- 
culated  from  starting  and  resulting  functions  shall  differ  less  than 
0.  Refinement  to  self-consistency  is  necessary  to  realize 


the  full  accuracy  of  which  Slater* s  method  may  be  capable.  However, 
the  functions  represent  the  first  phase  of  such  a  calculation*  More- 
oever,  since  they  differ  little  from  Jucys*  functions  (which  by  virtue 
of  arising  from  a  variational  principle  represent  the  "best"  one-electron 
approximation  embodying  the  Pauli  principle),  they  must  possess  consider¬ 
able  accuracy*  Therefore,  since  they  arise  from  a  plausible  unipotential, 
they  may  be  used  as  they  stand  in  calculations  requiring  unipotential* 

The  method  is  not  more  tedious  or  time-consuming  than  the  Hartree- 
Fock  method  (indeed,  it  is  less  so  by  reason  of  the  unipotential);  and 


7 


yields  a  more  accurate  set  of  wave  functions  than  the  Hartree  method, 
for  the  quadrivalent  carbon  stem  at  least*  Nevertheless#  carrying  the 
solution  through  to  self-consistency  is  a  difficult  project  for  the  lone 
worker;  especially  so,  because  the  self-consistent  procedure  is  not  nec- 


19 

esearily  convergent,  so  that  various  devices  must  be  used  to  obtain  con¬ 


vergence. 


19.  3.  Pratt,  Phys.  Rev*,  88,  1219,  *52. 
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APPENDIX  I 


An  expensive  literature  exists  on  the  Hartree-Fock  equations*^7  As 
4 

Modified  by  Dirac,  the  sf  equations  for  the  atonic  orbitals  read: 


6  J  1  \  ’  ’ 

To  put  this  in  the  form  of  a  Schrftdinger  equation.  Slater  rewrites  the 

exchange  term,  i«e«,  the  second  summation,  thus: 

(b)  1  1* rt  =  £ 


This  device  puts  (a)  in  the  Schoftdinger  form  with  an  effective  potential 
for  the  ^ith  electron  given  by 

(•)  -|+L  -  £ 

The  terms  of  this  potential  are  interpreted  physically  by  S  later  ^-^t/r 

is  simply  the  electrostatic  field  of  the  nucleus;  the  first  summation 

accounts  for  the  screening  effect  of  the  “charge  aloud"  of  electrons 

iurrounding  the  nucleus,  for  u*(2)u‘(2)  is  simply  the  SchrSdinger  charge 
•  ^  > 

density  for  the  Igth  electron,  and  the  integral  its  Coulomb  effect;  the 
second  summation  provides  "•••correction  due  to  the  fact  that  the  elec¬ 
tron  really  does  not  act  on  itself*"^  The  exchange  term  alone  of  the 
effective  potential  is  not  the  same  for  all  electrons;  to  render  it  so. 
Slater  strikes  a  mean  of  exchange  terms  weighted  with  the  probability  that 

the  electron  with  the  coordinates  (1)  is  in  the  state  _ii 

2: 

<(i)Uth)  / Z? 


<a Six)  <Jj( x)<Xf(^V 

dZx- 


17 •  J*  C.  Slater,  Phys«  Rev,,  V5.  210,  *50; 

V.  Fock,  Z*  Physlk,  £l,  12$7 

D*  and  W#  Hartree,  Proc>  Roy«  Soc«,  A 150  <  9*  '55» 

V.  also  footnotes  1,  5#  and  H#  ft®  well  as  references  given  by 
Slater,  loc •  primo  cit • 

18.  Slater,  Phys«  Rev>,  82*  *51« 


- 


c  v  '  - 


o 


4  *0 


.1 


s. 


1  0 

'1.  >?  0  I  J*  11  .:i'l  *<t 

' 

■il~  '  '  ...  •  '  ■  '•  :  ‘  ” 

•  :  t  _  i  ' 

■>.  .  . 

J  ■■  •  i  n  ~ 

/ 


.  ...  i 

i  t _  _  *  _  c  •  ’  *  ■ 

^  t  ?.  *  »  vJ  d 

t 

*  V  _  v*.  -  *  •*  * 

t  t  t-  ♦  _ 

*  y  _  * 

5  t  < _  <* _ -  ^  ‘ 


d 


— 24- 


The  result  i© 
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ror 


However,  in  view  of  the  insensitivity  of  the  core  electrons  to  the  pre¬ 
sence  of  other  electrons,  (ir*  sec.  1,  p*2),  a  better  approximation  to 
the  truth  obtains  if  the  valence  electrons  are  supposed  to  move  in  an 
exchange  which  results  from  averaging  (b)  over  them  alone.  With  the 
indices  as  assigned  in  sec*  2,  this  exchange  is 

The  denominator  works  out  to  jj>l(2s|r)  v  5PX(2p\rX]  /4  7Tr^  while  the 
summation  written  out  in  full  is 

(I,  ♦!  41  4-1  A I  )  4.  (I  +141*1  4 1  41  4 1  4-1  4T*r  421  421  ). 

N  v*  -iC  n'  V  «4  >f  »<•  4*-  TS  *>«*  At  4C  S* 

The  sums  in  parentheses  are  B(r)/4iv  and  0(r)/4rr  respectively,  essentially 
the  functions  tabulated  in  Table  I*  Combining  the  above  results  yields 

W(r)  *  [B(r)+C(r)]/(>l(2.|r)+5Pi(2p|r71  . 
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APPENDIX  II 


Evaluation  of  I  (y*) — 3y  definition 


=  ^  \  :>l(1)uI(2)*(1/rxl)u$(2)ufe(l) 


Now 

u:(l)u*(2)ua(2)u.(l) 

:  (5/K5ipP,'  P(lstr)p(8p»r)P(le\r)E(2pt^)/rfr‘, 

where  cos^r,  and  the  sinea  have  been  written  in  terms  of  the  associated 
Legendre  polynomials  ?J(co30). 

To  perform  the  integrations,  using  the  orthogonality  properties  of 
the  exponentials  and  Legendre  polynomials,  it  is  necessary  to  invoke 
the  following  expansions,  given  by  A*  Somraerfeld  in  Partial  Differential 
Equations  yi  Physics,  Pages  l^t  ff  • ;  (  ^c<x\enr^*c  New 

Vrj=  !Vr(;f%n(co?7  rjr*’  x  C  r, 

“  ( 1/r)  2  Pn ' c osry )  ( r,  /rj ,  r %  ^  r, 

P n(cos7/  =  2  r^rJ,Y  ‘^}e  ^  *■;  (n  -Iml  .  i.  (n  ♦  tn\)  l, 

where  cosy  -  1^,i^/r|r^_,  and  P  (£)  is  the  Legendre  Polynomial  of  degreen* 
The  volume  of  the  integration  must  be  divided  into  two  regions,  one  with 
raY  ri  *  and  the  °^^er  with  r%<  rx,  in 
converge  uniformly* 

The  angular  integration  for  the  region  r,_<  r%  reads: 


.  °<»  t  ( xr' 


'n-lml  !.'(n+  trt 

■Jjdr/p,'  ((Jp;  (£)  ( ^r)  ( rl/r)Epyf,)pl*’J^n  -l)l/(n*  m  )! 

r  ,  o**.f 

-  (2lvr,)  (r^/r  IP,  (fj  I  ,  1  "on  using  th  .ty  properties.  To 

average  this  over  angles,  it  was  (l/4*C)  X  and  integr 


ed  over  all  angles; 
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the  other  region,  the  roles  of  r  and  r_iaust  be  reversed. 

,,  now  bobmes  /<*>  ^ 

_  POMO  P(^Wi)  j_  Y  (is  o_b  V  ^ 

The  results  of  such  evaluationaStt  given  in  sec.  2  above. 
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